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ABSTRACT 

A Chandra ACIS-S imaging observation of the nearby galaxy M81 (NGC 3031) reveals 9 luminous 
soft X-ray sources. The local environments, X-ray spectral properties, and X-ray light curves of the 
sources are presented and discussed in the context of prevailing physical models for supersoft sources. It 
is shown that the sample falls within expectations based on population synthesis models taken from the 
literature though the high observed luminosities (L & s ~ 2 x 10 36 to ~ 3 x 10 38 erg-s -1 in the 0.2-2.0 keV 
band) and equivalent blackbody temperatures (T e ff ~ 40 to 80 eV) place the brightest detected M81 
objects at the high luminosity end of the class of supersoft sources defined by previous ROSAT and 
Einstein studies of nearby galaxies. This is interpreted as a natural consequence of the higher sensitivity 
of Chandra to hotter and more luminous systems. Most of the sources can be explained as canonical 
supersoft sources, accreting white dwarfs powered by steady surface nuclear burning, with X-ray spectra 
well-fit by hot white dwarf local thermodynamic equilibrium atmosphere models. An exceptionally bright 
source is scrutinized in greater detail as its estimated bolometric luminosity, L\, i ~ 1.5 x 10 39 erg s _1 , 
greatly exceeds theoretical estimates for supersoft sources. This source may be beyond the stability limit 
and undergoing a phase of mass outflow under extreme conditions. Alternatively, a model in which the 
observed X-ray spectrum arises from an accretion disk around a blackhole of mass -^OO/lcosz) 1 / 2 M Q 
(viewed at an inclination angle i) cannot be excluded. 

Subject headings: binaries: symbiotic — stars: evolution — stars: atmospheres — novae — white 
dwarfs — X-rays: stars 



nearby galaxy M81 reveals nine supersoft source candi- 
dates identified by broadband X-ray colors (§0). Four of 
these sources are located within the old population of bulge 
stars and all but one of the remainder are coincident with 
M81 spiral arms. A search of known objects finds two of 
the bulge sources are spatially coincident with sources of 
enhanced [OIII] A5007 A emission: consistent with neb- 
ulae photoionized by the X-ray sources. No other spa- 
tial correlations were discovered in searches of multiwave- 
length archival data, contemporary radio and Ha images, 
or published source catalogues. The two brightest super- 
soft source candidates are visible in archival ROSAT PSPC 
or HRI X-ray images. 

Detailed X-ray spectral analysis of the three brightest 
sources is presented in Their effective temperatures, 
based on line-blanketed local thermodynamic equilibrium 
white dwarf atmospheres and blackbody model fits, span 
T e //~50 to 80 eV. The corresponding inferred bolomet- 
ric luminosities are L& ;~3 x 10 38 to ~1.5 x 10 39 erg s _1 
ranking the brightest M81 objects among the hottest and 
brightest of all previously-known supersoft sources. Ex- 
trapolating the spectral fits to the fainter sources implies 
the observed luminosities of the M81 supersoft source pop- 
ulation range from L b s ^2 x 10 36 (near the sensitivity 
limit) to -3 x 10 38 erg-s" 1 in the 0.2-2.0 keV band. The 
brightest sources show signs of X-ray variability over the 
duration of the Chandra observation and over the approx- 
imately 10-year light curve observed by ROSAT (§|4|). 
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1. INTRODUCTION 

Luminous supersoft X-ray sources have effective black- 
body temperatures of 15-80 eV and bolometric luminosi- 
ties 10 36 to 10 38 erg s" 1 (Kahabka & van den Heuvel 1997). 
Lacking significant emission above ~0.5-l keV, supersoft 
sources are easily identified using even low-resolution X- 
ray detectors. However, only 10 Galactic objects are listed 
as supersoft sources in the latest compilation of Greiner 
(2000a) while population synthesis studies (Di Stefano & 
Rappaport 1994, Yungleson et al. 1996, Di Stefano & Nel- 
son 1996) predict of order 1000 supersoft sources could 
be active in ordinary galaxies like our own. Local objects 
are therefore of limited use for supersoft source population 
studies. The best opportunity is to explore nearby galax- 
ies at high galactic latitude where intervening absorption 
is low, all the sources are at a common distance, a single 
pointing can detect a large number of sources, and their 
spatial distribution can be measured. Indeed, surveys con- 
ducted using the ROSAT and Einstein observatories have 
identified 8 supersoft sources in the LMC, 4 in the SMC, 
and 34 in M31 (see Greiner 1996 for a detailed listing). The 
high sensitivity and moderate spectral resolution make the 
Chandra X-ray Observatory an ideal facility to systemat- 
ically extend supersoft source population studies to other 
nearby galaxies. This is the rationale for the present in- 
vestigation. 

A 50 ks Chandra ACIS-S imaging observation of the 
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The X-ray signatures of the M81 sources are compared 
to well-studied objects in the (rather heterogeneous) class 
of supersoft sources in an effort to identify possible physi- 
cal mechanisms responsible for the observed emission (§||). 
The brightest source cannot easily be explained by prevail- 
ing theory while the other two sources for which reliable 
spectral fits can be made show strong similarities to the hot 
close-binary accreting white dwarf (WD) systems CAL 87 
and CAL 83, respectively. Some prospects for future ob- 
servations are discussed (§||) and a comparison is made 
to the supersoft source populations observed in M31 with 
ROSAT (Supper et al. 1997, Kahabka 1999) and in M101 
by Chandra (Pence et al. 2001). 

2. OBSERVATION OF SUPERSOFT SOURCES IN M81 

A 49926 second observation of a portion of the galaxy 
M81 was obtained with the Chandra Advanced CCD 
Imaging Spectrometer (ACIS) on 2000 May 7. The nu- 
cleus and bulge of the galaxy were located near the center 
of the back- illuminated S3 device of the spectroscopy array 
operating in imaging mode. The observation was taken in 
faint timed exposure mode at 3.241 s-frame _1 . The global 
properties of the sources detected in this data are given 
in Tennant et al. (2001). Standard Chandra X-ray Center 
processing has applied aspect corrections and compensated 
for spacecraft dither. A charge transfer inefficiency (CTI) 
corrector algorithm (Townsley et al. 2000) 5 was then ap- 
plied to the Level 1 event list to partially correct for the 
charge loss and charge smearing effects of CTI in the ACIS 
detectors. The data were then cleaned of bad pixels and 
columns. We selected the standard grade set and events 
in pulse invariant (PI) channels cor responding to ~ 0.1 to 
8.0 keV for source detection (§ |t|) and -0.2-2.0 keV for 
spectral analysis (§ ||). The 50 ks exposure corresponds 
to a limiting supersoft source observed 0.2-2.0 keV lumi- 
nosity of — 2 x 10 36 ergs-s _1 for a 3er signal-to-noise ratio 
assuming the weak supersoft source candidates have the 
same spectrum as the brightest object in the class. 

The detector viewing area covers 57% of the optical ex- 
tent of the galaxy, defined as the ellipse of major diameter 
26/9 corresponding to the D 2 5 diameter as tabulated in 
de Vaucouleurs et al. (1991), oriented at position angle 
157°, and with major-to-minor axis ratio 2:1 correspond- 
ing to the 60° inclination angle of M81. This area includes 
all of the back-illuminated (BI) S3 device, approximately 
half of each of the front-illuminated (FI) S2 and S4 de- 
vices, and the outer corner of the FI device 13. The data 
from each device is analyzed independently due to differing 
energy resolutions, low-energy responses, and background 
signals. 

In addition to this primary dataset, a 2.4-ks ACIS-S 
image taken 2000 Mar 21 and numerous ROSAT PSPC 
and HRI datasets were used to construct long-term light 
curves of the brightest sources (see Immler & Wang, 2001 
for details of the ROSAT observations). 

2.1. Source Identification Criteria 

X-ray sources were identified using the source detec- 
tion method described in Tennant et al. (2001) that com- 
pares the data to an analytic function that approximately 
matches the telescope's point-spread function (PSF) and 

5 Matching response matrices were also provided by L. Townsley 



includes off-axis broadening. Small spatial regions cen- 
tered on each source with a size encompassing the 95% 
encircled energy radii of the PSF (at 1.5 keV) were then 
selected for analysis. These regions are < 2" diameter 
on-axis and exceed 10" only for sources ^8' off-axis. 

Background-subtracted source counts were binned into 
three broad bands defined as S (0.3-1.0 keV), M (1.0- 
2.0 keV), and H (2.0-8.0 keV) and the X-ray colors 
MS = (M - S)/(M + S) and HS = (H — S)/{H + S) 
were constructed. For this purpose, two background spec- 
tra were defined for sources on S3: One for the bulge where 
contributions from the wings of the PSF of the bright nu- 
cleus and from a relatively high level of unresolved X-ray 
emission was previously detected in the data (Tennant et 
al. 2001) and another for the disk. Separate background 
spectra are defined for sources on S2 and on S4. In all 
cases, regions enclosing all identified sources are first ex- 
cluded from the data then the spectrum of the remaining 
area was assigned to the background. Figure |l| shows the 
resulting background-corrected color-color diagram for all 
81 sources detected on S3 above a signal-to-noise ratio 
of 3.5 (Tennant et al. 2001). The bulge and disk back- 
ground spectra are also represented. Some faint sources 
far from the aimpoint are background-dominated. Their 
background-subtracted count rates in certain energy bands 
can be negative and the absolute values of their resulting 
colors may exceed unity, as seen in Figure |l|. 
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Fig. 1. — X-ray color-color diagram for M81 sources on the BI de- 
vice S3. Background-subtracted source counts are binned into three 
bands; S (0.3-1.0 keV), M (1.0-2.0 keV), and H (2.0-8.0 keV). 
Points labeled Disk and Bulge represent the colors of the background 
spectra in these two regions. 

Supersoft sources lie near (—1,-1) in this figure (no 
emission above 1.0 keV). Those sources with color values 
(MS < —0.5, HS < —0.5) were tenatively classified as su- 
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persoft sources. This criteria selects 10 sources on S3, two 
sources on S2, and none on the remaining CCDs. 

The list of sources was compared to catalogues of ob- 
jects to eliminate those with properties inconsistent with 
supersoft sources. Supernova remnant candidate no. 17 in 
the catalogue of Matonick & Fesen (1997) lies within the 
positional uncertainty of one of the supersoft source can- 
didates detected on S3. Though this source is weak, its 
spectrum appears more similar to the other X-ray sources 
associated with supernova remants than to the other su- 
persoft source candidates. Matonick & Fesen (1997) re- 
port an [OIII]/Ha ratio for the supernova remnant can- 
didate of 0.6 and an [SII]/Ha ratio of 1.2. Rappaport et 
al. (1994) predict far different ratios, 2.4 - 6.1 and 0.1 - 
0.5, respectively, for ionized regions surrounding super- 
soft sources. This source is therefore excluded from the 
list of candidate supersoft sources even though its colors 
(MS = -0.88 ± 0.39, HS = -1.10 ± 0.45) lie within our 
selection criteria. 

Both candidate sources on S2 are coincident with bright 
(mv^Q rnag) GO stars in the PPM catalogue of stars. 
While they have X-ray colors (MS, HS) ~ (-0.6, -0.8) 
they are excluded from further consideration. 

The positions of the remaining candidate supersoft 
sources were compared to optical features in Digitized 
Sky Survey (DSS) and archival Hubble/WFPG2 images 
to search for potential uncatalogued foreground objects. 
Late- type stars with active coronae and AM Her-type mag- 
netic cataclysmic variables (CVs) are the only foreground 
objects with X-ray colors falling within our selection cri- 
teria. Foreground late-type stars are easily identified in 
optical images as their X-ray flux is typically a small frac- 
tion of their optical flux. Am Her-type mCVs have a soft 
component (T ~ 10-100 eV) that may dominate any hard 
bremsstrahlung (kT ~ 10 keV) component and appear su- 
persoft. However, their K or M star companions should be 
visible optically out to distances well beyond the Galac- 
tic disk in the direction of M81. It remains possible that 
some of the weakest supersoft source candidates could be 
very bright CVs, with late M-type companions, and escape 
optical detection provided they are located some ~7 kpc 
above the plane of the Galaxy. 




The observed properties of the remaining 9 candidate 
supersoft sources, all located on the BI S3 device, are 
listed in Table 1 and their positions are shown superim- 
posed on the second generation DSS image of M81 in Fig- 
ure |[ Hereafter, individual supersoft sources in M81 will 
be referred to in order of observed brightness, column 1 of 
Table 1, as source Nl through source N9. 

Table 1 

Observed supersoft sources in M81 



No. 


RA 


DEC 


Flux" 


location 




(2000) 


(2000) 






1 


09 55 42.15 


69 03 36.2 


74.95 


bulge 


2 


09 56 08.96 


69 01 06.6 


10.24 


arm 


3 


09 55 53.00 


69 05 20.3 


3.30 


arm 


4 


09 55 28.38 


69 02 44.6 


1.09 


bulge 


5 


09 55 37.58 


69 03 16.2 


0.95 


bulge 


6 


09 55 55.97 


69 03 12.5 


0.78 


bulge 


7 


09 56 14.16 


69 02 26.2 


0.63 


disk 


8 


09 55 47.92 


68 59 28.2 


0.59 


arm 


9 


09 55 48.13 


68 59 15.5 


0.41 


arm 



"background-subtracted, 0.2-2.0 keV, 



Fig. 2. — Digitized Sky Survey image of M81 with the positions of 
the 9 X-ray-detected supersoft source candidates superposed. Image 
is ~15xl5 arcmin. Circles have 10" radii. 



in units of 10 cts-s 



The lack of candidate supersoft sources on the FI de- 
vices is attributable to the small fraction of M81 falling 
within the viewing field of the FI devices, to the scarcity 
of X-ray sources (of all types) found far from the central 
regions of M81, to the weaker low energy response of the 
FI devices, and to the decreased off-axis source sensitivity. 
The detection limit for supersoft sources on the FI devices 
is ~ 1.2 x 10 37 erg s _1 or about 6 times brighter than the 
corresponding limit on the BI device S3. 

2.2. Local Source Environments 

None of the candidate supersoft sources appear ex- 
tended in the Chandra data. As listed in Table 1, can- 
didate sources are found throughout the bulge and disk 
of M81 and several lie within ^200 pc of spiral arms 
(1' = 1.086 kpc for the adopted distance of 3.6 Mpc to 
M81, Freedman et al. 1994). 

While several of the candidate supersoft sources are in 
crowded, source-rich regions of this well-studied galaxy, 
only two can be unequivically identified with catalogued 
objects observed at other wavelengths based on positional 
coincidence. These are the bulge sources N4 and N5. 
Three of the 4 bulge supersoft source candidates lie within 
the -4' observing field of the [OIII] A5007 A survey of 
M81 planetary nebulae (PNe) performed by Jacoby et al. 
(1989). Source N4 is coincident with object ID 68 (ap- 
parent magnitude 25.33 at A5007 A) and source N5 is co- 
incident with object ID 116 (apparent magnitude 24.78) 
in the list of Jacoby et al. (1989). The observed X-ray 
luminosities of these sources exceeds that of typical plan- 
etary nebulae (10 30 to 10 32 erg s _1 ) by several orders of 
magnitude assuming they are at the distance of M81. The 
one known exception is 1E0056. 8-7154, a supersoft source 
associated with the Small Magellenic Cloud PNe N67 
(Wang 1991) with a bolometric luminosity 2 x 10 37 erg s^ 1 
(Heise, van Teeseling, & Kahabka 1994). The lack of 
detectable emission at other wavelengths, particularly in 
ffubbic/WFPC2 images, indicates they are not foreground 
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objects lying, by chance, along the line of sight. Instead, 
the [OIII] A5007 A emission could come from regions ion- 
ized by the X-ray source in a way similar to that observed 
in the LMC supersoft source CAL 83 (Remillard, Rap- 
paport, & Macri 1995). The flux observed by Jacoby et 
al. (1989) is consistent with the supersoft source ioniza- 
tion model predictions of R appap ort et al. (1994). This is 
discussed in more detail in §5.2.4. 



There arc no point-like objects observed at other wave- 
lengths spatially coincident with any of the remaining su- 
persoft source candidates in M81. It is noteworthy that 
the exceptionally X-ray bright source Nl is within the field 
observed by Jacoby et al. (1989), and the recent survey of 
Magini et al. 2001, but is not seen at [OIII] A5007 A. This 
source lies in a featureless region of the bulge dominated in 
optical bands by a high density of unresolved bulge stars. 
There is no nearby object discernable in archival Hub- 
We/WFPC2 images of the region, nor in UV images cen- 
tered at A2490 A and A1520 A (Hill et al. 1992), a recent 
(2001 May 21; A. Shafter, private communication) CCD 
image taken through a 70A wide Ha filter, the continuum- 
subtracted Ha data of Devereux et al. (1995), the A6 and 
A20 cm radio continuum observations reported by Kauf- 
man et al. (1996), or recent A6 cm data (T. Pannuti 2001, 
private communication) . 

3. X-RAY SPECTRAL ANALYSIS 

The spectra of all 9 candidate sources are shown in Fig- 
ure They typically show a rise from low energies to a 
peak at ~0.45 to ~0.7 keV followed by a rapid decline to- 
ward higher energies. The shape of the low-energy portion 
of the spectra is determined by a combination of extinction 
and a decreasing ACIS detection efficiency at low energies. 
Above the peak, the shape is determined by a roughly ex- 
ponential decline of the intrinsic source spectra. 

In spectra such as these, the inferred bolometric lumi- 
nosity is much higher than that observed with the peak of 
the intrinsic spectral energy distribution occuring at en- 
ergies lower than observed. The inferred luminosities are 
very sensitive to the model parameters and to calibration 
uncertainties in the low-energy response of the BI devices. 
Though there is often significant flux below 0.2 keV, data 
at energies below this value were excluded from the fitting 
because of these uncertainties. 

A sufficient number of counts for spectral analysis were 
accumulated from only the three brightest sources listed in 
Table 1. The spatial regions selected for spectral analysis 
were chosen centered on each source as described above 



(§2.1). Accompanying background spectra were extracted 
for each source from annular regions centered on the source 
with inner and outer radii of 30 and 50 pixels, respectively, 
and containing of order 200 counts. Note this procedure 
differs from the large background regions selected for con- 
structing the broadband colors described previously. The 
procedure used in this section more accurately represents 
local variations in background. Spectra have been grouped 
into spectral bins containing a minimum of 20 counts and 
fitted using the XSPEC spectral fitting package (Arnaud 
1996) to blackbody and to model atmosphere spectra of 
hot white dwarfs calculated in the local thermodynamic 
equilibrium (LTE) approximation. 

The X-ray emission from supersoft sources is believed 



to originate from nuclear burning on the surface of a WD 
star. The nuclear burning takes place at high optical depth 
within a geometrically thin, high density atmosphere with 
the resulting emergent spectrum highly modified by the 
cooler outer layers. Under these conditions, hot WD LTE 
model atmosphere spectra is a more appropriate approxi- 
mation than is a simple blackbody assumption. Therefore, 
we developed tables of plane-parallel LTE model atmo- 
spheres as an alternative model for spectral fitting. Under 
certain circumstances the WD atmosphere may become 
extended or even partially ejected as a result of rapid nu- 
clear burning. Spectra of extended, expanding, spheri- 
cal atmospheres tend to be flatter than that from plane- 
parallel hydrostatic atmospheres with the same color tem- 
perature at maximum flux. Model atmospheres including 
an optically thick wind are currently in preparation. 
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Fig. 3. — Spectra of nine supersoft source candidates. Labels de- 
note source number and count rates over the 0.2—2.0 keV Chandra 
energy band (Table 1). Data have been binned by 8 PI channels 
(~ 117 eV). Vertical axis units are counts/sec. 



The fundamental parameters of our plane-parallel LTE 
model atmospheres are the effective temperature, T e ff, 
surface gravity, g, and luminosity scale factor K = 
L/{AircrT^fD 2 ) where L is the bolometric luminosity, 
tr is the Stefan-Boltzmann constant, and D is the dis- 
tance from the source to the observer. The model atmo- 
spheres were constructed following standard temperature- 
correction procedures as outlined in Mihalas (1978) in 
order to satisfy the hydrostatic and radiative equilib- 
rium constraints. The 15 most abundant elements from 
H through Ni are included with bound-free opacities com- 
puted from the photoionization cross sections of Verner et 
al. (1993,1995). Line blanketing is included using ^1200 
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of the strongest spectral lines from the CHIANTI (v. 3.0) 
atomic database (Dere et al. 1997). 

A series of LTE model atmosphere spectra were con- 
structed in this fashion spanning the range of parame- 
ters 10 5 < T eff < 1.3 x 10 6 and 7.5 < \og(g) < 9.5, 
as appropriate for WD atmospheres, with the scale fac- 
tor K used as the model normalization in XSPEC. Tables 
of these models were used in conjunction with a photo- 
electric absorption component in fitting to the observed 
spectra. A separate table was constructed for each of two 
values of photospheric metal abundances, Z/Zq = 1 and 
Z/Zq — 0.01, relative to the solar composition given by 
Anders & Grevesse (1989). The lower metal abundance 
is roughly compatible with the observed M81 metallicities 
(~0.03) recently reported by Kong et al. (2000). In ad- 
dition, a He-rich model consisting of Z/Zq — 0.03 metals 
and 99% He by weight was also constructed. 

For reference, computed hot WD LTE model atmo- 
sphere spectra at several temperatures are shown in Fig- 
ure [|. Many of the edges seen in these spectra will prove 
key diagnostics of the properties of the candidate supersoft 
sources. 




100 200 500 1000 



Energy (eV) 

Fig. 4. — Hot white dwarf LTE atmosphere model spectra (solid) 
and corresponding blackbody spectra (dash) at four effective tem- 
peratures and a solar abundance composition. From top to bottom, 
the temperatures are T e ^/10 5 K = 10, 8, 6, and 4. The surface 
gravity is log(g)=9 except for the hottest model where log(g)=9.5 to 
avoid exceeding the Eddington limit luminosity. Prominent bound- 
free edges are apparent at 490 eV (C VI), 670 eV (N VII), 740 eV 
(O VII), and 870 eV (O VIII). 

3.1. Source Nl 

The spectral fit results for the brightest supersoft source 
candidate, source Nl, are listed in Table 2. This and sub- 
sequent tables for the other bright sources include one col- 
umn for each model. Rows list the \ 2 fit statistic, the 
fitted model parameters Nh (hydrogen column density), 
T e ff (blackbody or effective temperature), log(g) (surface 
gravity), E e d ge (absorption edge energy) and the accom- 
panying r (edge optical depth) . All models include photo- 
electric absorption from a solar abundance column. In the 
case of model zvarabs, the abundance of C and O in the 
absorbing column are allowed to vary. The table lists the 
enhancement of these elements relative to the solar values. 
Also listed are the observed luminosity, L b s , in the 0.2- 
2.0 kcV band, the unabsorbed luminosity, Lx, in the same 



band and the bolometric luminosity, Lb i, as scaled from 
the model normalizations assuming a distance of 3.6 Mpc 
to M81. 

The observed spectrum, best-fit absorbed blackbody 
model with added absorption edge at 0.26 keV (second 
model of Table 2), and fit residuals for source Nl are shown 
in Figure [| A simple blackbody model is statistically su- 
perior to either of the LTE WD atmosphere models be- 
cause of the lack of absorption edges in the observed spec- 
trum. However, there is a feature in the blackbody fit 
residuals near 0.3 keV. This feature can be modelled as an 
edge at 0.26 ± 0.01 keV (second column of Table 2) which 
is very close to the neutral C K-edge at 0.277 keV. Since 
there is no reason to expect the galactic absorbing column 
to have exactly solar abundances, we also tried a variable 
abundance model. If a small "redshift" (z=0.04) is ap- 
plied, this model also provided an excellent fit to the data 
as shown by the third column in Table 2. It is worth not- 
ing that the ACIS energy scale is uncertain at the lowest 
energies and thus the feature is likely a Carbon edge and 
the "redshift" is due to this uncertainty. As the ACIS re- 
sponse is also uncertain at these energies we cannot deter- 
mine whether the excess Carbon is intrinsic to the source 
in M81, lies within our galaxy, or is a detector effect. 
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Fig. 5. — Upper panel: Observed spectrum of source Nl (sym- 
bols) and best-fit blackbody model with solar-abundance photoelec- 
tric absorption and added absorption edge at 0.26 keV (solid line). 
Lower panel: Fit residuals. 

LTE models with an additional Raymond-Smith ther- 
mal emission component are an improvement over the sin- 
gle LTE component models because this added thermal 
component tends to fill the ~0. 7-1.0 keV region of the 
model spectrum underestimated due to the presence of 
absorption edges in the LTE component. While the phys- 
ical basis for such a thermal component has precedence in 
modeling of, e.g., the recurrent nova U Sco (Kahabka et al. 
1999) and the symbiotic nova SMC 3 (Jordan et al. 1996), 
and there is reason to believe source Nl may be near a state 
of dynamical outflow (see discussion in §[|), the fact that 
a simple blackbody model provides an equally-acceptable 
fit allows one to conclude that the thermal component is 
not necessary. It should be noted that ^2% of the incident 
flux at 0.5 keV is piled up in the spectrum of source Nl 
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Table 2 
Model Fits for Source Nl 



Model 


blackbody 


blackbody 


blackbody 


LTE 


LTE 


LTE + ray 


diskbb 






+ edge 


w/ zvarabs 6 


Z = Z e 


2 = O.O1Z 


z = z Q 




Xi/daf 


1.40/47 


0.98/45 


1.05/44 


2.75/46 


3.55/47 


1.24/44 


1.61/47 


Nh/W 20 cm -2 


r 7+0.7 
°-'-0.6 


A O+i-O 




4 9+0.6 

^• z -0.4 


a n+ - 6 

*- u -0.4 


4-4±g;| 


7 4+0.8 
'•^-0.8 


T e// eV 




01 +3 
8II4 




^n+1 

69IJ 


n 4 +2 


66I{ 


woti 


-E ed9e keV, r 




0.26l°° 2 , 1.6 












C,0 enhancement 






6.611:?, i.5§:! 










log(g) cm 2 s _1 






7+0.04 
°- ' -0.04 


o+0. 05 
S -°-0.03 


8.5i°;f 




thermal fcT, /(Lx) c 












0.35^° 07 12 , 0.18 




Lobs/W 38 erg s" 1 


n n +0.02 
°- u -0.02 


n +o.o4 

J - u -0.04 


^• y -0.03 


9 q+0.02 
z ' 3 -0.02 


9 n+0.04 
z - a -0.04 


q n+0- 01 

°- u -o.oi 


q n +0. 03 
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''variable-abundance absorption column with best-fit redshift z = 0.04 
c Raymond-Smith model temperature (keV) and fraction of Lx due to this component 



resulting in detected events near 1 keV. The observed flux 
at 1.0-1.1 keV is 3.5% of that in the 0.5-0.55 keV band. 
Thus some 60% of the observed flux at 1.0-1.1 keV is due 
to pileup. 

He-rich LTE model atmospheres were also fitted to this 
spectrum but no statistically-acceptable fits were obtained 
and they are omitted from Table 2 for brevity. The larger 
He opacity near the Hell edge in these models leads to a 
slightly lower temperature and a slightly steeper temper- 
ature gradient in spectrum-forming (O.l^Tij^Sl, where tr 
is the Rosseland mean opacity) regions of the atmosphere. 
The emergent flux therefore decreases near the Hell edge 
(due to lower temperatures) and increases slightly at the 
peak of the spectral energy distribution (due to the steeper 
gradient). The overall effect in the Chandra observable en- 
ergy band is to produce a steeper spectrum than the H-rich 
LTE models. The spectrum is further steepened by the 
bound-free edges from intermediate-mass metals resulting 
in a poorer model of the data. 

The blackbody shape of the observed spectrum of source 
Nl is reminiscent of that expected from a geometrically 
thin optically thick accretion disk surrounding a black hole 
(Shakura & Sunyaev 1973). For completeness, therefore, 
an accretion disk emission model was also fitted to the 
observations. This model accumulates the spectrum from 
an optically thick accretion disk observed at inclination 
angle i from a superposition of blackbody spectra from 
concentric annular elements representing the radial depen- 
dence of the disk temperature profile (Mitsuda et al. 1984, 
Makishima et al. 1986). The model parameters (model 
diskbb, Table 2) are the normalization, proportional to the 
innermost disk radius, Ri n , for a known inclination an- 
gle and distance to the source, and the temperature at 
Ri n . Assuming Ri n corresponds to the last stable Keple- 
rian orbit (Makishima et al. 2000) implies a mass for the 
central object. The best-fit normalization for source Nl 
corresponds to a black hole mass Ms^^l200/(cosi) _1 / 2 
assuming the last stable orbit is at 3 Schwarzschild radii. 
The best-fit temperature at Ri n , Ti„^100 eV, implies a 
bolometric luminosity Lfc o ;~10 39 erg s _1 . Not surpris- 
ingly, the color temperature Ti n is slightly higher than 
that obtained from the blackbody models. The shape of 
the blackbody disk model is a power law of photon index 



—2/3 at low photon energies steepening to the Wien por- 
tion of a blackbody of temperature T\,b = 0.7T;„ at higher 
energies (Makishima et al. 1986). Thus the combination 
of an innermost disk radius temperature slightly higher 
than the blackbody model temperature and a higher ab- 
sorption column density is just what is needed to make the 
disk blackbody model mimic a standard blackbody spec- 
trum. This is especially true in view of the fact that only 
the steep Wien portion of the spectrum lies within the 
Chandra energy band for low values of Ti n and the power 
law portion of the blackbody disk model spectrum lies en- 
tirely below the low-energy instrumental cutoff. The lack 
of detectable emission at longer wavelengths at the posi- 
tion of source Nl (§|]^) and the lack of an observable hard 
power-law tail argues against the blackbody disk model. 

The luminosity of source Nl is exceptionally high. The 
observed 0.2-2.0 keV luminosity, L b s ~3 x 10 38 erg s _1 , 
is slightly above the Eddington limit for a 1.4 M H- 
accreting star. The unabsorbed luminosity on the same 
energy range and the inferred bolometric luminosities are 
much higher though with high uncertainties. In contrast, 
L bs is tightly constrained and its value does not vary 
among the models (Table 2). Another source of uncer- 
tainty, independent of model fit uncertainties, is the gain 
behavior in the BI device at low photon energies and at the 
-120°C focal plane temperature in use at the time of obser- 
vation. Very nearly 1/3 of the observed flux from source 
Nl falls in the 0.2-0.4 keV energy band. A 20% uncer- 
tainty in the flux below 0.4 keV corresponds, therefore, 
to an additional 7% uncertainty in L b s - Conservatively, 
then, the estimated bolometric luminosity of source Nl is 
at least several times 10 38 erg s _1 and, depending on the 
model, may be as high as a few 10 39 erg s _1 . This result 
places severe constraints on physical models for this source 
(§§)•, 

Tt is notable that the large range of Lx and L^oi is not 
obviously reflected in a large range of the model fit pa- 
rameters Njj and T e ff. From Table 2, the range of these 
parameters among the various models (excluding diskbb) 
are N H = 4.9 ± 2.0 x 10 20 cm" 2 and T eff = 75 ± 11 eV. 
This substantiates the previous conjecture that the in- 
ferred luminosities are very sensitive to small variations in 
the model parameters. However, the differences in shape of 
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the blackbody and LTE model atmosphere spectra do not 
give systematically large differences in the derived tem- 
peratures, column densities, and luminosities as is often 
reported for supersoft sources. 

The hot WD LTE model parameter \og(g) is directly 
related to the mass of the assumed white dwarf. The re- 
sulting white dwarf mass is estimated to be 0.9 to 1.1 M Q 
for source Nl. 

3.2. Source N2 

The spectral fit results for source N2 are listed in Ta- 
ble 3. The observed spectrum, best-fit LTE atmosphere 
model, and fit residuals are shown in Figure ^. 

Table 3 

Blackbody Model Fits for Source N2 
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blackbody 
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19.8 +8 7 °, 
98±if 
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the presence of highly-ionized O in a very hot atmosphere. 
The LTE models provide significant improvements over 
the blackbody model. The results are not very sensitive to 
the assumed metal abundances though the 1% solar metal 
abundance model provides an improved fit. At the derived 
temperatures, T eff ~ 67-69 eV, the hot WD LTE mod- 
els produce numerous absorption edges. When convolved 
with the instrument response, the most prominent edges 
are those from highly ionized O at 740 eV (O VII) and 870 
eV (O VIII). In the low-abundance model, the edges are 
not as deep and an additional Ne IX edge (1196 eV) ap- 
pears in the model but is statistically insignificant in the 
fit. 

The observed 0.2-2.0 keV luminosity derived from the 
hot WD LTE model atmosphere spectrum is much less 
than the Eddington limit for a hydrogen-accreting WD. 
The inferred range of bolometric luminosities exceed this 
limit by a factor of 2 to 3. The white dwarf mass inferred 
from the LTE models is 0.7 to 0.9 M Q . 

3.3. Source N3 

The spectral fit results for source N3 are listed in Ta- 
ble 4. The observed spectrum, best-fit hot WD LTE model 
atmosphere spectrum, and fit residuals are shown in Fig- 
ure 0. 

Table 4 
Model Fits for Source N3 



Model" 



blackbody Z = Z© Z = 0.01Z G 



X'i/dof 
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10.7+^i 6 
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Fig. 6. — Upper panel: Observed spectrum of source N2 (symbols), 
best-fit LTE WD model atmosphere spectrum with 1% solar abun- 
dance metals and solar-abundance photoelectric absorption (solid 
line). Lower panel: Fit residuals. 

In contrast to source Nl, a blackbody model is clearly a 
poor fit to the data. A blackbody with an added absorp- 
tion edge at 0.85±°;°3 keV gives a better fit suggesting 




0.3 0.5 1 

channel energy (keV) 

Fig. 7. — Upper panel: Observed spectrum of source N3 (symbols), 
best-fit WD LTE atmosphere model spectrum with solar abundance 
metals and solar-abundance photoelectric absorption (solid line). 
Spectrum and model have been rebinned to show structure. Lower 
panel: Fit residuals. 
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Unlike the stronger sources, the spectrum of source N3 is 
fit equally well by all three models and there is no spectral 
feature to meaningfully discriminate between models. In 
particular, the lack of flux above ^0.7 keV means that the 
dominant absorption edges present in the LTE atmosphere 
models are the relatively weak edges at C VI (490 eV) and 
N VII (670 eV). In this case, a slightly more absorbed 
blackbody model relative to the LTE models adequately 
mimics the observed steep decline above ~0.5 keV. 

The observed 0.2-2.0 keV luminosity, L a b s , is less than 
the Eddington limit for a hydrogen-accreting WD and Lboi 
is comparable to the Eddington limit though, again, it is 
highly uncertain. The white dwarf mass inferred from the 
WD atmosphere LTE models is ^0.5 M Q for the solar 
abundance model and is effectively unconstrained by the 
1% solar metallicity model. 

3.4. Other Sources 

Statistically constrained model fits to the remaining 
sources could not be achieved. The observed count rates 
and similarities to the spectra of the brighter sources al- 
low a crude extrapolation of the derived luminosities of 
the brighter sources to the remaining sources by the fol- 
lowing reasoning: The ratio Lb i/L bs ~ 5 for source NI 
and ~ 10-20 for sources N2 and N3. The latter two sources 
are associated with spiral arms. They have higher absorb- 
ing columns than source Nl which accounts for the higher 
Lboi/L b s ratio. An estimate of Lboi can be made by as- 
suming the ratio Lb i/L b s = 15 for sources on the spiral 
arms and = 5 otherwise and noting that L D b s is propor- 
tional to the observed count rate (Table I). Applying this 
approximation to the remaining sources using the data 
from Table 1 gives the following values of Lboi m units of 
10 37 erg s" 1 : 2.2 (N4), 1.9 (N5), 1.6 (N6), 1.3 (N7), 3.6 
(N8), and 2.4 (N9). The uncertainties in these estimates 
are at least a factor of three since the actual absorbing 
columns are unknown. This uncertainty is, alas, compa- 
rable to some of the uncertainties in Lboi obtained for the 
brightest sources. 

Alternatively, luminosities can be estimated from the 
count rate and assuming a (blackbody) temperature and 
a column density. The weakest supersoft source detected 
has an observed luminosity of 2 x 10 36 erg s^ 1 assuming 
an 80 eV blackbody spectrum and a column of Nh = 8x 
10 20 cm -2 . This corresponds to Lboi ~ 1.3 x 10 37 erg s _1 . 
Lboi increases to 5.2 x 10 37 erg s _1 if the assumed temper- 
ature is lowered to 50 eV. 

Temperatures cannot be so easily estimated but they 
are most probably of order 40-50 eV as higher tempera- 
ture sources are expected to be rare based on theoretical 
models of supersoft sources (see §H) and cooler sources are 
difficult to detect with Cha.ndra-A.CIS. Even an intrinsi- 
cally bright {Lboi = 10 38 erg s _1 ) blackbody source with 
T e ff = 25 eV and moderate absorbing column (2 times 
Galactic) produces only 24 counts in the BI device and 
would be undetectable in the FI chips. 

4. X-RAY TIMING ANALYSIS 

4.1. Short-Term Variability 

The 50 ks Chandra light curves of the brightest 3 super- 
soft source candidates are presented in Figure pi Detected 
events have been binned on 1000 sec intervals. 




TIME (sec of 2000-5-7) 

Fig. 8. — Chandra light curve of the brightest supersoft source 
candidates binned on 1000 sec intervals. All standard-grade events 
in the 0.2-2.0 keV bandpass, within valid good-time intervals, and 
within the same spatial regions as used for spectral analysis, are 
included. Error bars represent lcr statistical uncertainties. Back- 
ground contributions are negligible: Typical background count rates 
are 10 — 6 counts s — 1 pixel -1 . Labels denote source number and av- 
erage count rate obtained by fitting a constant to the data (shown 
as horizontal lines through the data). 



The observation was of sufficient duration that the 
Kolomogorov-Smirnov statistic can be evaluated for the 
brightest sources to test the hypothesis that the sources 
are constant. The test shows sources Nl and N2 are vari- 
able at high confidence but the remaining sources contain 
insufficient data to be conclusive. Power spectra were also 
generated for the brightest sources to search for pulsations 
or other periodic behavior. Only source Nl had sufficient 
signal to obtain significant results. No periodicity was de- 
tected. 

An additional 2.4 ks ACIS-S imaging observation of M81 
obtained on 2000 Mar 21, 47 days previous to the 50 ks ob- 
servation, was acquired through the Chandra X-ray Cen- 
ter data archive. Source Nl is the only supersoft source 
candidate with a sufficient number of counts in this short 
exposure for analysis. The light curve of source Nl during 
this observation, binned on 300 sec intervals, is shown in 
Figure ^. Also shown are the corresponding light curves in 
the energy bands 0.2-0.5 and 0.5-1.0 keV and the time de- 
pendence of the background during the observation. The 
source brightness declines from a maximum at the begin- 
ning of the observation to a statistically insignificant level 
within ^900 s and remains at this low level for the remain- 
der of the observation. No other (bright) source displayed 
similar changes in X-ray flux during the 2.4 ks observation. 
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Fig. 9. — Chandra 2ks light curve of source Nl binned on 300 sec 
intervals. The uppermost curve includes all standard grade events 
in the 0.2—2.0 keV band, the next lower curve represents the 0.5- 
1.0 keV band (dashed) followed by the 0.2-0.5 keV band (dotted). 
The lower curve represents the nearby background scaled by the ra- 
tio of the source to the background area. Error bars represent la 
statistical uncertainties. 

The dramatic drop in flux suggests an eclipse event. The 
observation interval is too short to establish a brightening 
of the source signifying the end of the eclipse (the observa- 
tion ended 100 sec into the last interval and this interval 
contains 4 counts, one of which is probably a background 
photon based on its PI value). No similar feature is de- 
tected in the longer Chandra observation and no similar 
feature can be seen in the relatively poor signal-to-noise 
ROSAT data discussed below. 

Assuming the drop in flux signals the onset of an eclipse 
of an ~1 M Q compact object, that Roche-lobe overflow 
from the companion is occuring, that the eclipse duration 
exceeds 1600 s, and that no eclipse occurs during the 50 ks 
exposure, weak constraints on the orbital elements of any 
binary system can be derived. These constraints are easily 
met by massive companions (^3 M©) that have evolved off 
the main sequence. Systems in this mass range have peri- 
ods greater than ~1.6 days and eclipse fractions of order 
15%. Missing an eclipse is then quite likely though wit- 
nessing the onset of an eclipse in a short exposure would be 
fortuitous. Furthermore, such com panion stars should be 
visible in the Hubble images (§|2.2|), A more likely scenario 
for supersoft sources is a slightly-evolved main sequence 
donor in the mass range 1.3-2.5 M© (van den Heuvel et 
al. 1992) or a "helium Algol" with a few-solar-mass com- 
panion at the onset of He accretion (Iben & Tutukov 1994, 
Yungelson et al. 1996). Orbital periods range from 10 hr to 
30 hr for the main sequence donor scenario and <20 hr for 
the He-accreting systems (Yungelson et al. 1996). There 
is at best a one-in-three chance that no eclipse is observed 
during a 50 ks (13.9 hr) observation of systems with a 30 hr 
period and the probability decreases for the shorter-period 
systems. Thus, though an eclipse cannot be excluded by 
the X-ray observations it is not very likely. 

Assuming the source is exiting an eclipse phase during 
the last ~300 s of the short exposure, then the eclipse lasts 
only of order 1800 s. The orbital period must then be of 



order a few hours only since ~15% of the period is spent 
in eclipse. A period this short could not be missed in the 
50 ks observation. 

Alternatively, nuclear burning on the surface of a WD 
may be occuring through weak shell flashes that either 
cause expansion which quenches burning or expels mate- 
rial which obscures the underlying atmosphere as it cools 
with more violent flashes followed by a larger flux decline. 
This may explain the Figure || light curve where the X-ray 
flux is roughly a factor-of-two brighter during the onset 
of the observation than the average value observed during 
the 50 ks observation (Figure ||). Fitting an exponential 
to the first 1000 s of this light curve gives an e-folding 
timcscalc of ~ 350 s. This is comparable to the sound 
crossing time over the surface of a WD for reasonable es- 
timates of the sound speed in the burning layer. A change 
in accretion rate may also trigger large excursions in the 
X-ray brightness. Southwell et al. (1996) have suggested 
this mechanism may explain the factor of ^20 rise in X- 
ray flux observed (Schaeidt, Hasinger, & Trimmer 1993) 
in the supersoft source RXJ 0513.9-6951. However, the 
timescale for contraction of the photosphere leading to 
the required change in T e ff is of order days (Livio 1992). 
RXJ 0513.9-6951 is probably a low-mass WD system with 
a correspondingly low X-ray luminosity (Greiner 2000b) 
far below that observed for source Nl. 
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Fig. 10. — Chandra 50 ks light curve of source Nl in the 0.2- 
0.5 keV (top) and the 0.5-1.0 keV (middle) energy bands binned on 
1000 sec intervals. The bottom panel displays the ratio of the two 
energy bands. Error bars represent 1<t statistical uncertainties. 

Interestingly, the spectrum is substantially harder and 
the brightness is higher during the bright portion of the 
2.4 ks observation compared to the 50 ks observation. 
There is even some flux above 1.0 keV in the 2.4 ks 
data. The mean photon energy during this bright phase 
is ~730 eV while during the 50 ks observation the mean 
is ^475 eV (Figure ||). Part of the spectral hardening, 
including most of the flux above ~1.1 keV, is caused by 
pileup in the high rate data. Nevertheless, the spectrum 
is definitely harder (or, hotter in the context of blackbody 
spectral models) during this bright phase. There is no 
clear signature of increasing absorption as the source en- 
ters the low rate region of the 2.4 ks light curve as would be 
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expected in the onset of an eclipse. No other source shows 
a similar harder spectrum during the 2.4 ks observation. 

The other alternative discussed previously in the context 
of spectral models for source Nl is X-ray emission from an 
accretion disk around a mid-mass black hole. Variability 
on many timescales is observed in both massive black holes 
(AGNs) and in stellar-mass black holes (X-ray transients). 
Variability is generally caused by instabilities in the disk 
structure and changes in the accretion flow and is usually 
accompanied by spectral evolution among high soft states 
and low hard states. The 50 ks light curve of source Nl 
in each of two energy bands is shown in Figure [l^ along 
with the ratio of the soft band to the hard band. This 
shows that source Nl varies differently in different energy 
bands but it does not help distinguish between standard 
supersoft source and mid-mass black hole scenarios. 

4.2. Long-Term Variability 

The fluxes observed by Chandra suggest the brightest 
4 or 5 supersoft source candidates may be detectable in 
archival ROSAT observations of the field. M81 was ob- 
served 9 times by ROSAT/PSPG over the interval 1991 
March through 1994 April and 11 times by ROSAT/HRI 
from 1992 Oct through 1998 April. Source Nl is present 
in ROSAT/HRI images but is too near the bright nucleus 
of M81 to be resolved in ROSAT/PSPC data. Source N2 
is clearly present and isolated from other X-ray sources 
in both PSPC and HRI observations. Source N3 is con- 
fused with a bright nearby object in both ROSAT instru- 
ments. Source N4 is likewise too near the nucleus for pos- 
itive detection. Source N5 is marginally detected at the 
2.8<7 level in the longest of the HRI exposures. This cor- 
responds to a Chandra luminosity of ~ 2 x 10 37 erg s _1 ; 
consistent with a constant source luminosity between the 
two observations. The remaining sources are too weak 
to be detected at a reasonable confidence level. The en- 
tire ROSAT dataset has been independently analyzed by 
Immler & Wang (2001). Chandra source Nl is within 3. "2 
of HRI source H25 of Immler & Wang (2001) and is listed 
as a variable source with no other comment. Source N2 
is within 3. "3 of HRI source H36 which is coincident with 
PSPC source P44. Source P44 is a soft source based on its 
PSPC hardness ratios (Immler & Wang 2001). No other 
Chandra-detected supersoft source candidate is among the 
ROSAT sources identified by Immler & Wang (2001). 

The combined ROSAT/HRI, PSPC, and Chandra light 
curves of sources Nl and N2 are displayed in Figure [tl]. 
Here, the average count rate of each individual observation 
is shown combined to form each data point. The Chandra 
and ROSAT/PSPC observed count rates have been scaled 
to ROSAT/HRI observed count rates using the PIMMS 
tool (Mukai 1993) and assuming spectral properties for 
the two sources as derived in section ||. This scaling can 
account for some of the disparity seen, for example, in 
the lower panel of Figure |l| (source N2) where the PSPC 
data and the two Chandra points are lower than the in- 
tervening ROSAT/HRI fluxes. Nevertheless, both sources 
are clearly present in the data spanning >7.5 yr. As dis- 
cussed in the next section, this fact limits the possibility 
that these sources could be classical novae. 
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Fig. 11. — Combined ROSAT and Chandra light curves for Source 
Nl (upper panel) and N2 (lower panel). Count rates have been 
scaled to ROSAT/HRI count rates assuming the spectral param- 
eters derived in §0 from the 50 ks Chandra data. ROSAT/HRI 
data are marked by squares, ROSAT/PSPC data with triangles, 
and Chandra data with circles. Errors denote lcr statistical uncer- 
tainties. Labels denote average observed count rates. 

5. DISCUSSION 

Supersoft sources are a well-established class of X-ray 
emitting objects. Of the 19 well-studied systems (e.g., 
Greiner 2000a), all within the Galaxy or the Magellenic 
Clouds, nine are accreting WDs in close binaries (CBSS) 
that are nuclear burning accreting material in a steady 
state, four are classical or recurrent novae, three are sym- 
biotic systems, and one is a planetary nebula nucleus. The 
common source of X-ray emission among these sources is 
nuclear burning cither of accreting material (in a steady 
state or in a nova flash) or of residual fuel (following, e.g., 
the formation of a planetary nebula) . 

The assumption that all supersoft sources result from 
nuclear burning on WD stars is adopted here for purposes 
of discussion. Two questions to be addressed are: (1) Is 
the observed population and its distribution within M81 
consistent with this theory and (2) Are the observed X- 
ray properties of the individual sources consistent with 
previously-observed members of this rather heterogenous 
class (and, if so, which ones)? 

5.1. The M81 Population of Supersoft Sources 

Nine supersoft source candidates have been identified in 
M81 based on broad-band X-ray colors. All are located 
on the S3 device. Four of these are located in the bulge of 
M81 and the remainder are aligned with the spiral arms 
with the exception of one object in the inter-arm disk re- 
gion. The bulge comprises 30% of the area of the S3 de- 
vice 6 thus 2.7 of the 9 supersoft sources should occur in 
the bulge if the spatial distribution were uniform on S3. 
However, 54 of the 97 X-ray sources of all types on S3 are 
located in the bulge (Tennant et al. 2001) and therefore 5 
supersoft sources would be expected there if the supersoft 
sources follow the general distribution of X-ray sources. 
Twenty-one of the 43 disk X-ray sources are aligned with 
the spiral arms while 4 of the 5 supersoft disk sources are 
on the arms. 



6 The bulge is defined here as the region interior to the inner Lindblad resonance at radius 4±0.2 kpc (Kaufman et al. 1989 and references 
therein). This is larger than the bulge radius (2.55 kpc) adopted by Tennant et al. (2001) resulting in 13 more X-ray sources within the bulge. 
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Fig. 12. — Distribution of supersoft source systems over effective temperature (left) and bolometric luminosity (right). Shown are the 
predicted and observed distributions for the galaxy M81 (see text) and the distribution of known sources in the Galaxy and Magellenic Clouds 
as catalogued by Greiner (2000a). The M81 predicted distributions assume the sources are distributed normally over both temperature and 
luminosity. This normal distribution is the curve in the left panel. For clarity, the corresponding distribution is omitted from the right panel 
as it appears as a narrow asymmetric curve centered on 3xl0 37 erg s _1 in the logarithmic scaling used here. 



The observed properties of the supersoft sources in M81 
favor high mass WD systems (M]yn~1.0 M©) and hence 
a younger population of progenitor stars (zero-age main 
sequence mass ~7 M ). This favors locations near star 
forming regions such as spiral arms. Conversely, the in- 
terstellar absorption is relatively higher in the spiral arms 
than in the disk or bulge which would tend to obscure the 
lowest-temperature supersoft sources. Indeed, the two su- 
persoft sources located on the spiral arms with sufficient 
counts for spectral fitting (N2 and N3) have relatively high 
X-ray absorbing columns and temperatures of 50-60 eV. 

Objects in the bulge are expected to be an older pop- 
ulation (>9 Gyr, e.g., Kong et al. 2000) except in the 
central ~50" where an enhancement in recent star for- 
mation activity is suggested by UV (Hill et al. 1992, Re- 
ichen et al. 1994) and Ha (Devereux, Jacoby & Ciardullo 
1995) imaging and consistent with the concept of gas- 
infall toward the center of M81. This may have begun 
as recently as ^400 Myr ago as a result of an encounter 
with NGC 3077 (Thomasson & Donncr 1993) or M82 (do 
Grijs, O'Connell, Gallagher 2001). All the bulge supersoft 
source candidates lie outside this central region and appear 
to have relatively low absorbing columns consistent with 
the lack of gas and dust. The standard supersoft source 
model (van den Heuvel et al. 1992) requires companion 
stars in the mass range 1.3-2.5 M Q to maintain steady 
nuclear burning of accreting H on the surface of the WD. 
Stars of this mass range evolve from the main sequence on 
timescales of order several 10 s to greater than 10 9 years 
and so the required accretion could be occuring in the cur- 
rent epoch in the bulge if the donor stars were formed in 
an encounter. 

How many supersoft sources are expected in the 50-ks 
Chandra observation? Di Stefano & Rappaport (1994) es- 
timated the number of supersoft sources in external galax- 
ies observable by ROSAT by carefully accounting for the 
effects of interstellar absorption on an assumed parent pop- 
ulation taken from the work of Rappaport, Di Stefano, & 
Smith (1994). The parent population was that of systems 



whose luminosities and temperatures matched the steady- 
burning regime computed by Iben (1982). This omits 
several classes of supersoft source including wind-driven 
symbiotic systems, novae, and recurrent novae. Neverthe- 
less, it is possible to apply the basic approach of Di Ste- 
fano & Rappaport (1994) to the Chandra observation of 
M81. Scaling from their estimate of the total population 
of supersoft sources in M31 by the ratio of blue luminosi- 
ties [i.e., assuming supersoft sources contain massive WDs 
and hence evolved from a young stellar population (Motch, 
Hasinger, & Pietsch 1994, Di Stefano & Rappaport 1994)], 
a total of 1700 supersoft sources are expected in M81. 

To estimate the fraction detectable by Chandra, the 
population distributions in temperature and bolometric 
luminosity as derived by Rappaport, Di Stefano, & Smith 
(1994) can be roughly approximated by normal distribu- 
tions with means T = 24 eV and L = 3 x 10 37 erg s _1 and 
standard deviations of 13 eV and 8.5 x 10 36 erg s _1 , respec- 
tively (see also Fig. 3 of Di Stefano & Rappaport 1994). By 
randomly sampling values of T and L from these two in- 
dependent distributions, assuming a blackbody spectrum 
and a modest column of JY^ = 8 x 10 20 cm~ 2 (twice 
the Galactic value), and using the PIMMS software tool 
(Mukai 1993), the predicted Chandra observed count rate 
distribution can be accumulated. Assuming 20 counts in a 
50-ks observation as a detection threshold, only 1% of the 
supersoft sources in M81 should be detectable with the BI 
device S3. 

As the viewing field of S3 encompasses 23% of the D25 
area of M81, ~4 of the 1700 supersoft sources in the par- 
ent population should be detected in the 50-ks observation 
on S3. This is consistent with the observed population of 
9 sources in light of the large uncertainties involved in- 
cluding a factor-of-two uncertainty in the normalization of 
the parent population (Di Stefano & Rappaport 1994) and 
the neglect of other classes of supersoft source systems by 
these authors (Yungleson et al. 1996). 

The simulation also shows that only systems with tem- 
peratures T>40 eV and bolometric luminosities ~ 2 x 



12 



10 37 erg s _1 are likely to be detected by Chandra. This 
result can be compared to the distribution of temperatures 
and luminosities reported in Greiner (2000a) for 21 super- 
soft sources (two of which have not been optically clas- 
sified) and to the 9 supersoft sources discovered in M81. 
Figure 11 shows this comparison. (For brevity, it has been 
assumed here that the 6 weakest supersoft sources in M81 
all have temperatures of 45 eV, see §^J.) The distributions 
of the Greiner (2000a) sample are clearly broader than 
predicted or observed for M81. This is a consequence of 
the greater diversity of local environments in which these 
sources are found and the dissimilar instruments used to 
investigate them. The peak of the observed M81 distri- 
butions in temperature and luminosity coincide with the 
predicted distribution but these peaks are dominated by 
the weak sources and are, therefore, highly uncertain. The 
bright observed sources are, however, located at higher 
temperatures and luminosities than expected. 

5.2. The Nature of Individual M81 Supersoft Sources 
5.2.1. Source Nl 

By far the brightest supersoft source candidate in the 
sample lies within the bulge of M81 approximately 52" 
from the nucleus. Its spectrum is well-fit with a simple 
blackbody model of 78 - 86 eV (T ~ 9 x 10 5 K) with a 
modest absorption (Nh^6 X 10 20 cm~ 2 ) resulting in an 
implied bolometric luminosity exceeding 10 39 erg s _1 . Al- 
though the bolometric luminosity is highly uncertain, even 
the observed luminosity exceeds the Eddington limit for a 
Chandrasekar-mass WD. 

The brightest well-studied supersoft sources are typi- 
cally novae or recurrent novae with peak bolometric lu- 
minosities often exceeding 10 38 erg s -1 . One of the best 
examples is U Sco. Kahabka et al. (1999) investigated 
BeppoSAX observations of this recurrent nova taken ~20 
days after the peak of the most recent optical outburst. 
They find a temperature of ~ 9 x 10 5 K and a bolometric 
luminosity of up to 2 x 10 38 erg s _1 from a combination of 
non-LTE atmosphere models and an optically-thin ther- 
mal emission component, assumed to arise from a wind, 
added to fit the ~l-2 keV portion of the observed spec- 
trum. 

Although source Nl has X-ray properties similar to 
those of U Sco, it is unlikely that source Nl is a nova. 
Novae occur in systems with relatively low mass transfer 
rates where a layer of fuel is accumulated for a period of 
time before compressional heating causes ignition. The 
mass of the accumulated layer increases with decreasing 
WD mass. Consequently, the X-ray light curves evolve 
slowly, lasting up to 10 years or so, in systems with low- 
mass WDs. The luminosity is proportional to the WD 
mass M WD , L bo i~l.8 x 10 38 (M WD — 0.26) erg s™ 1 assum- 
ing H accretion, where Mwd is measured in solar units 
(Ibcn & Tutukov 1989). Thus the luminosity is low for 
these long-lived low Mwd novae and the recurrence time 
between events is very long. Conversely, systems contain- 
ing high mass WDs, such as inferred for U Sco, have a 
higher peak luminosity but are short-lived X-ray sources 
(~0.1 yr, e.g., Kato 1997). The recurrence times are also 
shorter but are still of order years and therefore inconsis- 
tent with the observed light curve of source Nl. (Note 
the short- and long-term variations observed in the X-ray 



light curve of source Nl are better described as transience 
or variability in contrast to recurrence as applied to no- 
vae.) 

At higher accretion rates is a regime of steady-state 
burning of accreting material. The maximum luminosity 
in the steady-burning regime is L^oi = 2.3 x 10 38 (Mwd — 
0.5) erg s _1 and occurs at the maximum accretion rate 
M= 7.8 x 10- 7 (M WD - 0.5) Mo-yr^ 1 (Iben 1982, again 
assuming H-rich accretion). The highest temperatures, up 
to ^85 eV, are obtained by the highest mass WDs as are 
the highest luminosities. However, an accretion rate of 
~5x 10~ 6 Moyr -1 is needed to obtain the bolometric lu- 
minosity inferred for source Nl assuming steady burning. 
At such high accretion rates, however, the photospheric 
radius expands to red giant dimensions; only a fraction 
of the fuel is consumed; and the X-ray emission, formed 
in a thin layer at the base of the envelope, is hidden by 
the overlying material. The burning may instead drive a 
strong, opaque wind (Hachisu, Kato & Nomoto 1996). If 
the outcome of a high mass transfer rate is formation of an 
envelope then this system should be observable optically 
as a hot giant or red giant star, contrary to the Hubble 
observations (§ |2.2| ) and not as a strong X-ray source. If 
the accretion drives a wind, then the photospheric tem- 
perature is only T ~ 10 5 K. However, if the accretion rate 
decreases the wind mass-loss rate will decrease as will the 
effective photospheric radius. Regions closer to the burn- 
ing layer will become exposed and the X-ray temperature 
and X-ray luminosity will increase dramatically (Hachisu, 
Kato, & Nomoto 1996). The subsequent evolution of sys- 
tems of this type were not followed by Hachisu, Kato, & 
Nomoto (1996). 

If the accreting material is He-rich, then the system 
can sustain a higher mass transfer rate while remaining 
in the steady nuclear-burning regime (Iben & Tutukov 
1989). The increase in mass transfer rate by roughly an 
order of magnitude translates into only a factor of 2 or 
so higher luminosity because of the lower specific energy 
generation rate of He relative to H so that the maximum 
luminosity is L ho ;^4.6 x 10 38 (M WD - 0.58) erg s _1 (Iben 
& Tutukov 1989). Although He-rich LTE white dwarf at- 
mosphere models failed to provide a satisfactory fit to the 
observed spectrum, He-rich acccretion remains an attrac- 
tive scenario for source Nl because the Eddington limit 
luminosity for He accretion is twice that for H. Helium ac- 
cretion can occur during a second phase of mass transfer to 
the compact star when the companion has exhausted He 
in its core (Iben & Tutukov 1994). The donor star in this 
scenario is either the nucleus of an early asymptotic giant 
branch or an evolved He-star remnant with a CO core. Un- 
der appropriate conditions, this mass transfer can be con- 
servative (Iben & Livio 1993) instead of lost in a common- 
envelope process (Iben & Tutukov 1985). The initial mass 
of the donor star tends to be rather large in this scenario, 
~6.5-9.5 M Q for steady burning, so that such systems 
should be rather rare in the old population of bulge stars 
in M81. Helium tends to burn explosively (Sion & Star- 
rfield 1993) unless the helium layer can be maintained at 
a high temperature (Jose et al.1993). A series of weak He 
shell flashes may occur for systems near the stability limit 
producing luminosities of order ^4x 10 38 erg s -1 , and tem- 
peratures ~70 eV (Iben & Tutukov 1994). These systems 
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tend to have less than ^20 hr orbital periods (Yungelson 
et al. 1996), comparable to or less than those of subgiant 
donors in H-accreting systems and of main sequence stars 
losing mass on a thermal timescale (~lday orbital peri- 
ods, van den Heuvel et al. 1992). Another channel for 
He accretion is a semidetached double degenerate system 
in which the secondary is a helium WD (Iben & Tutukov 
1989). Such a system must have an orbital period of order 
minutes (Verbunt & van den Heuvel 1995). 

In short, none of the standard nuclear-burning WD sce- 
narios predict luminosities as high as that inferred for 
source Nl. However, theoretical treatments have not rigor- 
ously explored the behavior of these systems in regimes be- 
yond the stability limit. Other alternatives exist. Neutron 
star binaries often exhibit a luminous thermal component 
though rarely as soft as observed here (and hence one of the 
original motivations for the nuclear-burning WD sceanrio 
for supersoft sources, van den Heuvel et al. 1992). One 
exception is the Be/X-ray binary RX J0059. 2-7138 with 
T e //~36 eV (Hughes 1994) though this source is accom- 
panied by a strong power law component extending the 
spectrum beyond the ROSAT band. Weakly magnetized 
neutron stars such as the X-ray burst sources, tend to have 
correspondingly weak power law spectral components but 
also have high characteristic temperatures during the burst 
phase (e.g., Lewin, van Paradijs, & Taam 1995). Similarly, 
typical black hole binaries display a "soft" component with 
characteristic temperature >1 keV and a hard power law 
component. 

A completely different scenario that may explain the 
observed X-ray properties of source Nl is X-ray emission 
from an optically-thick accretion disk surrounding a mid- 
mass black hole. The most attractive feature of this model 
is that the implied bolometric luminosities derived here 
are a small fraction of the Eddington limit. The observed 
X-ray spectrum is consistent with that from the inner- 
most disk radius assuming the central object has a mass 
M B h ~1200/(cosz) 1/2 . While no objects in this mass range 
have been confirmed, invoking the Eddington limit to X- 
ray-bright sources has led to suggestions of masses in ex- 
cess of 10-100 Mq, in many instances {e.g., Makishima et 
al. 2000) and as high as ~700 M© in the case of a lumi- 
nous source in M82 (Kaaret et al. 2001, Matsumoto et al. 
2001). 

5.2.2. Source N2 

Source N2 has a distinct edge at ^870 eV well fit by 
O VIII absorption in the LTE atmosphere model. This 
feature is commonly found in the hotter supersoft sources 
such as the LMC source CAL 87. Ebisawa et al. (2001) 
find a temperature of 75 eV for CAL87 from ASCA ob- 
servations with an absorbing column comparable to that 
obtained here for source N2 (Nh ~2x 10 21 cm -2 ). 

Parmar et al. (1997) report BeppoSAX observations of 
CAL 87 are equally well-fit with either blackbody, LTE, 
or NLTE atmosphere models with temperatures ranging 
from 42 eV for the blackbody fits to 57 and 75 eV for 
their LTE and NLTE models. The blackbody model im- 
plies L bol ~4 x 10 38 erg s -1 but the LTE and NLTE fits 
are consistent with more modest luminosities, ~3-5 x 10 36 
erg s . Including the O VIII absorption edge increases 
the blackbody temperature and decreases the bolometric 
luminosity to values similar to their NLTE models. 



This trend is opposite to the results reported here for 
source N2: blackbody models predict higher temperatures 
than the LTE models and bolometric luminosities compa- 
rable to or lower than the LTE models. The best-fit model, 
LTE with 1% of solar metal abundance, has a temperature 
of 67 eV and a hydrogen column density of 1.6 x 10 21 cm~ 2 
very similar to CAL 87. The observed luminosity is 
3.9 x 10 37 erg s _1 and L bo i~5.?>±^l x 10 38 erg s -1 , much 
higher than the models of Parmar et al. (1997) predict but 
consistent (within the large uncertainties) with the highest 
luminosity reported by Ebisawa et al. (2001) for CAL 87. 

CAL 87 is a close binary system that exhibits eclipses 
with a period of 10.6 hours. The eclipses are visible as 
shallow dips in the X-ray light curve. Small-amplitude 
variability of this kind cannot be distinguished in the rel- 
atively low-quality light curve of source N2. 

5.2.3. Source N3 

This source is also relatively highly absorbed but con- 
siderably cooler than the two brightest sources with 
T e ff ~ 50 eV. This value is more typi cal of many of the 
well-studied supersoft sources (Figure 5.1, Greiner 1996, 
2000a), such as the CBSS sources CAL 83 and IE 0035.4- 
7230. Parmer et al. (1998) report acceptable fits can be 
achieved with blackbody, LTE, or NLTE atmosphere mod- 
els applied to BeppoSAX observations of CAL 83. Their 
resulting temperatures are ~45 eV for the blackbody and 
LTE models and -33 eV for their NLTE fits. As with 
source N3 (§ |3.3| ), Parmer et al. (1998) report the low tem- 
perature and consequent lack of strong absorption edges 
does not allow a distinction to be made among the models 
(see also Figure |l|). BeppoSAX and ROSAT observations 
of IE 0035.4-7230 (Kahabka, Parmer, & Hartmann 1999) 
again do not allow either blackbody, LTE, or NLTE mod- 
els to be discounted though in this case a feature ascribed 
to C V and C VI absorption is evident. Temperatures de- 
rived by Kahabka, Parmer, & Hartmann (1999) range from 
-40±13eV from their blackbody fits to ~28±4 eV from 
their NLTE model fits to IE 0035.4-7230. Both CAL 83 
and IE 0035.4-7230 derived bolometric luminosities are of 
order a few times 10 37 , roughly an order of magnitude 
lower than the best-fit values derived here for source N3 
though within the large uncertainties in the latter. 

5.2.4. Source N4 and N5 

Though little can be said about the X-ray properties of 
these two sources beyond their broad-band classification as 
supersoft source candidates, these objects may be similar 
to the LMC source CAL 83. 

Source N4 and N5 are spatially coincident with objects 
ID 68 and 116, respectively, in the list of [OIII] A5007 A 
sources observed by Jacoby et al. (1989). This emission 
line is also prominent in the optical images of CAL 83 
(Remillard et al.1995) and is attributed there to ioniza- 
tion in the interstellar medium by a large UV and soft 
X-ray photon flux from the central source. Rappaport et 
al. (1994) modeled the ionized regions surrounding super- 
soft sources and calculated resulting optical line intensi- 
ties. The strongest predicted lines are the [OIII] A5007 A 
and Hell A4686 A lines. Adopting the foreground extinc- 
tion to M81 from Jacoby et al. (1989), their observed [OIII] 
A5007 A fluxes are consistent with those expected from the 
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supersoft source models. For example, the A5007 A flux 
from object ID 116 is ^4 x 10 -16 erg cm -2 s _1 while Rap- 
paport et al. (1994) predict 1.6 X 10~ 16 erg cm -2 s _1 and 
2.3 x 10~ 15 erg cm~ 2 s _1 from models with intrinsic lumi- 
nosities of 10 37 and 10 38 erg s _1 , respectively, and source 
temperatures of 4 x 10 5 K. 

Parmar et al. (1998) find a temperature T e //~30-50 eV 
for CAL 83 and L& ;^4 x 10 37 erg s _1 . These values are 
within the estimated values for sources N4 and N5. 

A PNe cannot be excluded based on the X-ray and [OIII] 
fluxes of sources N4 and N5. The central stars of PNe 
can have luminosities as high as that inferred for sources 
N4 and N5 but only the most massive central stars are 
hot enough to be detected in soft X-rays (e.g., Paczyn- 
ski 1971). This is the case for the SMC supersoft source 
1E00056.8-7154. This source is coincident with PNe N67 
(Wang 1991) but is anomalously X-ray bright for a PNe 
[T e// ~38 eV L bo i~2 x 10 37 erg s" 1 (Brown et al. 1994, 
Heise et al. 1994)]. Unfortunately, no obvious spectro- 
scopic distinction between PNe and supersoft sources has 
been identified (Di Stefano, Paerels, & Rappaport 1995). 

5.2.5. Other Sources 

The X-ray evidence in support of a supersoft source na- 
ture of the remaining candidates is based entirely on their 
X-ray colors. Their spectra and light curves are of insuffi- 
cient quality for detailed analysis. 

As there is no long-term X-ray light curve for these 
sources, their identification as novae cannot be ruled out. 
The rate of nova occurrence in M81 is rougly 20 yr -1 
(Delia Valle & Livio 1994) and, as discussed above, X- 
ray lifetimes range from months for the brightest sources 
up to ~10 years. For the values of L{, ; estimated in §|[ 
10 37 to 4 x 10 37 erg s , the corresponding WD masses are 
low, Mwd ~ 0.3 to 0.5 M Q . These masses correspond to 
long-lived novae, typically radiating as a (very soft) X-ray 
source for ~7 yr for M WD =0.6 M (Kato 1997). Visual 
magnitudes of novae can be as high as my ~ 22 to 23 
at the distance of M81 but the visible light rapidly de- 
clines after maximum light as the spectrum hardens and 
the X-ray flux rises. A 0.6 M Q WD nova, for instance, 
only becomes X-ray visible ^4 yr after optical outburst. 

The weak supersoft source candidates could also be 
steady-burning CBSS systems. The estimated L^oi cor- 
responds to a range of WD masses to <1.0. 
However, WDs in this mass range should be consider- 
ably cooler (T e //< 40 eV) than the 40-50 eV tempera- 
tures suggested by the instrumental selection effect. Be- 
low MvyD^0.65 M Q , the photospheric radius expands in 
the steady nuclear burning case and T e f f plummets (Iben 
1982). Thus only a narrow range of masses are possible 
for the weak supersoft sources if they are steadily burning 
their accreted hydrogen. 

Another evolutionary channel leading to the formation 
of a supersoft source is wind-driven accretion in symbi- 
otic systems (Sion & Starrfield 1994). These systems 
can undergo H shell flashes, reach luminosities of order 
~10 37 erg s _1 , and remain at this plateau for ^250 yr in 
some circumstatances. For many of the low-mass WDs 
modeled by Sion & Starrfield (1994), the flashes are weak, 
the envelope radius remains compact, and the resulting 
photospheric temperature is relatively high. However, the 
maximum temperatures achieved were T e //^20 eV. The 



supersoft source SMC 3 (RX J0048.4-7332) is a symbiotic 
nova of this type. Jordan, et al. (1996) derive a luminos- 
ity ib /^4 x 10 37 erg s _1 , a temperature of 22 eV, and 
a mass of ^0.8 M Q for this object. The weak supersoft 
source candidates are unlikely to be symbiotic systems be- 
cause of the difficulty detecting objects of such low tem- 
peratures in M81 with Chandra. An archival Chandra ob- 
servation of AG Dra, for instance, detects very little flux 
above ~ 0.3 keV. 

6. SUMMARY & PROSPECTS 
6.1. Chandra Observations of M81 

An important part of the ROSAT legacy has been the 
identification of supersoft sources as a phenomenologically 
distinct class of X-ray emitting object. The Chandra X- 
ray Observatory affords an excellent opportunity to extend 
the study of supersoft sources to other nearby galaxies. 

As a first step, the population of supersoft source candi- 
dates discovered in M81 has been investigated. The num- 
ber of supersoft sources detected and their spatial distri- 
bution are consistent with population synthesis estimates. 
X-ray spectral analysis shows the M81 sources are qual- 
itatively similar to previously-studied supersoft sources 
though the brightest M81 sources tend to be hotter and 
more luminous than is typical of the class. This was shown 
to be, in part, a natural consequence of the higher sensi- 
tivity of Chandra to the high-temperature region of the 
supersoft source distribution. 

Bolometric luminosities deduced here from the spectral 
fits often exceed the Eddington limit for spherical accre- 
tion onto a 1.4 M Q star. Such extreme luminosities lie 
beyond the previously-observed distribution of supersoft 
sources and, more importantly, cannot be reconciled with 
prevailing theories for supersoft sources as WDs powered 
by stable, steady-state, surface nuclear burning. But, the 
study of supersoft sources with the Chandra Observatory 
is still in its infancy. Important uncertainties remain. As 
in the ROSAT era, spectral modeling remains problematic. 
Simple blackbody models, like those used in this work, are 
known to greatly overestimate bolometric luminosities in 
some cases. LTE model atmospheres, also applied here, 
are appropriate when a well-defined photosphere is present 
such as in the high density gradient at the surface of a WD. 
Even this model fails when the nuclear burning is near the 
stability limit, the atmosphere becomes extended, and a 
wind drives mass loss as may be the case in the brightest 
supersoft source candidate in M81. In addition, the en- 
ergy calibration of the BI devices at low photon energies 
is uncertain and could also have a significant effect on the 
results reported here. 

Many of these uncertainties may be resolved or at least 
mitigated as pending Chandra observations of well-studied 
supersoft sources are analyzed and compared to previ- 
ous observations and to specific theoretical expectations. 
High-resolution grating spectroscopy will be indispensible 
in this regard as many of the calibration uncertainties can 
be circumvented while, at the same time, models for the X- 
ray spectra can be more tightly constrained. An example 
of the richness of such a spectrum is that of CAL 83 ob- 
tained by XMM-Newton/RGS (Paerels, et al. 2001). The 
population of supersoft source candidates discovered in 
M81 also reveals the value of isolating sources to inves- 
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tigate their local environments and of obtaining precise 
locations for future observations at X-ray and other wave- 
lengths. 

In spite of the current uncertainties, the remarkable 
properties deduced for the brightest supersoft source can- 
didate in M81 should not be overlooked. This source ex- 
hibits large variations in X-ray brightness on timescales of 
order hours yet maintains a high average X-ray brightness 
throughout the roughly 8 year span of combined ROSAT 
and Chandra observations of the region. In the context 
of the nuclear burning WD scenario, its high tempera- 
ture and luminosity imply a high mass, near the Chan- 
drasekhar limit for WDs, and a high accretion rate, above 
that for steady-state nuclear burning. While optical mon- 
itoring would be difficult, at best, future observations are 
strongly encouraged to more fully understand this source. 
The accurate location provided in Table 1 should aid sub- 
stantially in this regard. 

6.2. Observations of Other Galaxies 

The most extensive study of supersoft sources in an ex- 
ternal spiral galaxy similar to our own and M81 is the 
ROSAT survey of M31 (Supper et al. 1997, Kahabka 1999). 
Based on hardness ratios, Supper et al.(1997) identified 15 
supersoft source candidates not associated with SNRs or 
foreground objects in a 6.3 square degree M31 field. Ka- 
habka (1999), using a slightly modified selection criteria, 
identified an additional 26 candidates in the field. Ka- 
habka (1999) found a total of 7 of the combined 41 su- 
persoft sources are within the ^6 kpc bulge of M31 with 
the remainder evenly distributed over ~12 to 25 kpc radii. 
The lack of detected sources in the region ^6 to 12 kpc was 
explained as a consequence of a higher H column within 
this annulus. Kahabka (1999) suggests that the spatial 
distribution favors a disk population of younger stars by 
comparing to the populations of Cepheids and (older) blue 
stars. The supersoft source distribution does not follow 
that of novae which are bulge-dominated and follow the 
old stellar population found in the bulge. After accounting 
for possible foreground objects and SNRs, Kahabka (1999) 
finds 1 bulge source for every 4-7 disk sources. This is 
marginally less than reported here for M81 where 4 bulge 
sources were found for an extrapolated total population of 
9/0.57~16 supersoft sources within the D25 area of M81. 
Kahabka (1999) also estimates the temperatures for all the 
sources and finds they range from ^30-60 eV, typically, up 
to ^73 eV for the recurrent transient RX J0045. 4+4154. 
This latter source is also one of the brightest with an in- 
ferred i 6o( -10 38 erg s- 1 (White et al. 1995). 

A more detailed comparison between the supersoft 
source populations of M31 and M81 must await further 
analysis of the extensive Chandra M31 dataset. Only 
then can many of the instrumental selection effects be 
eliminated. In the interim, comparison to other Chandra- 
observed galaxies must suffice. 

Numerous deep observations of nearby galaxies have 
been performed by Chandra in the ^2 years since its 
launch. Discovery (or lack thereof) of supersoft source 
candidates have been reported in the literature for only a 
few galaxies. Sarazin, Irwin, & Bregman (2000,2001) re- 
port 3 supersoft sources in a sample of 90 point sources in 
the elliptical galaxy NGC 4697 to a limiting luminosity of 



~ 5 x 10 erg s based on hardness ratios. No super- 
soft sources were discovered in M84, an elliptical galaxy 
in the core of the Virgo cluster, based on the hardness ra- 
tio criteria (Finoguenov, Jones, & Kudritzki 2001). Pence 
et al. (2001) report the discovery of ten supersoft sources 
among 110 objects in M101 using color definitions differ- 
ent than those defined here. Pence et al. (2001) combined 
the spectra of the three most luminous supersoft sources, 
reportedly having similar spectral shapes, and find a best 
fit blackbody temperature of 72±2 eV and a mean (un- 
absorbed) luminosity of 1.4 x 10 38 erg s _1 . These values 
are similar to the brightest supersoft sources in M81 with 
the exception of source Nl which has an unabsorbed flux 
5-8 times higher. The remaining 7 supersoft sources in 
M101 have a distinctly softer spectrum. When combined, 
the best fit blackbody temperature is 47±2 eV and unab- 
sorbed luminosity is Lx = 1.1 x 10 37 erg s _1 for a column 
fixed to the Galactic value (N H = 1.2x 10 20 cm -2 ). Again, 
these values are consistent with those estimated here for 
the weak sources in M81 although, as shown above, a tem- 
perature of order 40-50 eV is expected on the basis of in- 
strumental selection effects and the shape of the parent 
population estimated from population synthesis models. 
The position of only one supersoft source candidate lies 
within the bulge of M101. However, M101 has a later 
Hubble type, Scd, and hence a smaller bulge than M81 
(Sab) extending only to ~0.'75, or approximately 3% of 
the S3 viewing field analyzed by Pence et al. (2001). 

These comparisons are intriguing. They suggest that 
perhaps there is another supersoft source population of 
anomolously X-ray bright and hot objects in nearby galax- 
ies similar to our own. Care must be exercised, of course, 
as there are, to date, only a handful of objects with these 
properties, no corroborating evidence in support of this 
conjecture from other wavebands, and no satisfatory the- 
ory predicting such extreme X-ray behavior. These com- 
parisons also suggest a trend across the Hubble sequence 
with fewer supersoft sources in early type galaxies com- 
pared to late types. The supersoft sources appear to be as- 
sociated with the younger population of stars found on spi- 
ral arms though some are found in the relatively old popu- 
lation of bulge stars. Irregular galaxies would then have a 
disproportionate number of supersoft sources. Di Stefano 
& Rappaport (1994) estimate 125 and 25 supersoft sources 
are active in the LMC and SMC, respectively, based on 
their blue luminosities relative to M31. Eight supersoft 
sources have been observed in the LMC and 4 in the SMC. 
This suggests either the population estimate is low or the 
observed fraction is high for these two galaxies. Of course, 
there is a strong selection bias that must be taken into 
account. 
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